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INTRODUCTION 


When  an  electromagnetic  wave  interacts  with  a  conductive  sphere,  it  induces  an 
eddy  current  which  is  detected  as  an  induced  diamagnetism.  Equations  for  the  mag¬ 
netic  polarizability  of  a  solid  sphere  have  long  been  available  in  the  literature1  and  in 
various  texts  on  electromagnetism,  for  example,  Landau  and  Lifschitz’s  Electrodynamics 
of  Continuous  Media.2  From  the  polarizability,  one  can  develop  expressions  for  the 
magnetic  permeability  of  a  spherical  particle  and  for  a  composite  of  such  particles  in 
an  insulating  binder.  While  difficulties  are  evident  when  the  calculated  composite  per¬ 
meabilities  are  compared  with  experimental  data,  the  formulas  are  at  least  useful  guides 
to  the  anticipated  behavior  of  such  composites.3 

Metal  coated,  insulating  spherical  shells  will  show  an  analogous  diamagnetism  in 
their  interactions  with  electromagnetic  waves.  However,  comparable  equations  for  the 
polarizability  and  permeability  in  the  literature  have  not  been  found,  although  they  may 
exist.  In  this  report,  it  shall  be  shown  how  such  formulas  may  be  derived  following  the 
approach  in  Landau  and  Lifschitz,  and  the  dependence  of  the  permeability  on  coating 
thickness  and  conductivity,  as  well  as  on  the  diameter  of  the  sphere  will  be  explored. 


THE  PERMEABILITY  OF  A  SPHERICAL  SHELL 


Assumptions 

A  few  assumptions  are  necessary  for  the  formulation.  Outside  the  shell,  it  is 
assumed  that  the  wavelength  is  much  greater  than  the  physical  dimensions  of  the 
particle  and  that  we  thus  have  a  quasistatic  field  without  sources  X  3  =  0). 

Inside  the  shell  (of  conductivity  a).  Maxwell’s  equations  and  the  appropriate  wave 
equations  for  the  magnetic  field  and  current  density  (  T  )>  assuming  a  sinusoidal  varia¬ 
tion  in  fi,  are  as  follows: 


?xE  = 

C 

(1) 

$•§  =  o  =  ^-5 

(2) 

^  ^  j*l  _  4jt T  _  4 no& 

(3) 

?x?xl  -  4*i,'"‘7  J 

<T 

(4) 

V2  ft  +  k2  ft  -  0, 

(5) 

1.  FORD,  G.  W„  and  WERNER,  S.  A.  Helicon  Oscillations  in  a  Sphere.  Physical  Review,  v.  8,  1973,  p.  3702-3709. 

2.  LANDAU,  L.  D.,  and  LIFSCHITZ,  E.  M.  Electrodynamics  of  Continuous  Media.  1st  ed.,  Addison  -  Wesley  Publishing 
Company,  Inc.,  Reading,  MA,  p.  193-194,  1960.  Also  see  2nd  ed„  Pergamon  Press,  New  York,  1984,  p.  205-206. 

3.  RUSSELL,  N.  E.,  GARLAND,  J.  G.  and  TANNER,  D.  B.  Absorption  of  Far  Infrtred  Radiation  by  Random  Metal  Ptrticle 

Composites.  Physical  Review,  v.  32,  1981,  p.  632-639. 


1 


where  k2  =  —  —  ^  °  ,  or  k  =  j-  where  <5,  the  skin  depth,  is  defined  as  3  =  x  a>  a ' 

0 

It  is  noted  that  we  have  assumed  that  Ohm’s  law  is  valid  in  the  static  limit;  i.e., 
the  relaxation  time  for  the  electrons  is  much  shorter  than  the  inverse  of  the  frequency. 
This  is  true  for  frequencies  below  the  mid-infrared  region.  Also,  the  local  assumption 
about  Ampere’s  law  presupposes  that  the  mean  free  path  of  the  electrons  within  the 
conductor  is  small. 

Formulation  and  Solution  of  the  Boundary  Value  Problem 

Consider  a  spherical  shell  of  inner  radius  b  and  outer  radius  a,  with  a  unit  exter¬ 
nal  magnetic  field  incident,  the  direction  of  which  can  be  taken  to  be  along  the  z  axis 
without  loss  of  generality.  Inside  the  inner  conductor  (r  <  b),  the  field  must  be  well 
behaved  at  the  origin,  and  must  be  a  (frequency  dependent)  constant;  therefore: 


A 

z  , 


(6) 


where  y  is  the  internal  field  strength. 


Outside  the  shell,  we  expect  a  dipolar  term  to  be  added  in  opposition  to  the  inci¬ 
dent  field  and  this  induced  dipolar  field  must  vanish  at  infinity;  therefore: 


ft  =  ~  [3  cos  9?  -  z  ]  +  z  ,  (7) 

where  a  is  the  magnetic  polarizability. 

Since  the  external  field  is  uniform  along  the  z  axis,  it  can  only  induce  dipole  terms 
within  the  shell,  and  J  will  have  only  an  azimuthal  component.  Solutions  to  Equation  4 
for  the  current  density  must  therefore  be  of  the  form 


sin  kr 


cos  kr 


]  z  X  r  =  f(r) 


/V 

Z 


/V 

r 


(8) 


from  which  we  can  determine  H  within  the  shell, 
ft,  =  (y  +  l^f^z  -  k2^  cos0r\ 


(9) 


where  f'  denotes  the  derivative  of  f(r)  with  respect  to  r. 

Applying  the  boundary  condition  that  the  field  be  continuous  at  r  *  a  and  r  =  b, 
the  following  equations  are  obtained: 
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at  r  =  b, 


3  y)  -l-  k2  f(b)  =  0, 

+  t2  f(b)  =  y, 
and  at  r  =  a, 

3  f  M  +  k2  f  (a)  =  > 

a  a3 

and 

+  k2  f  (a)  =  +  1 . 

a  a3 


(10) 

(11) 


(12) 


(13) 


For  convenience,  we  introduce  the  following  definitions:  B  =  kb,  A  =  ka, 

D  =  a3,  G  =  y  b3,  St  =  sinkb,  S2  =  sinka,  Cj  =  coskb,  C2  =  coska,  Tx  =  tankb, 

T2  =  tanka,  and  T  =  tan(A-B).  With  these  definitions,  Equations  10-13  become: 

P\  [3  (-  Si  +  B  CO  +  B2  Si]  +  p2  [3  (Ci  +  B  SO  -  B2  Q]  =  0  (14) 

Pi  [-  Si  +  B  Ci  +  B2  SO  +  p2  [Ci  +  B  Si  -  B2  CO  =  G  (15) 

pi  [3  (-  S2  +  A  C2)  +  A2  S2]  +  pj  [3  (C2  +  A  S2)  -  A2  C2]  =  -  3  a  (16) 

Pi  [-  S2  +  A  Q  +  A2  Si]  +  p2  [Cz  +  A  S2  -  A2  CJ  =  -  a  +  D.  (17) 

With  some  relatively  straightforward  algebra  (eliminate  Pi  and  p2  between  the 
equations  and  collect  terms),  one  obtains  solutions  for  the  polarizability, 

3  1  T  [f3  -  A2 )  (3  -  B2 )  +  9  A  B]  -  3(A-B)  (3  +  AB) 

&r  a2  T  (3  -  B2  )  +  3B 

and  for  the  internal  field, 


Y  = 


_ 3A _ 

(3-B2 )  sin  (A-B)  +  3B  cos  (A-B) 


(19) 


It  can  be  readily  verified  that  the  result  for  the  polarizability  reduces  to  the 
correct  formula  for  a  solid  sphere. 

The  field  outside  the  sphere  is  that  of  a  uniformly  magnetized  sphere  of  permeability 
/ip,  where 


H  = 


1  + 


1  - 


8jio 

~r 

Ana 

~1~ 


(20) 
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Calculated  values  of  /Up  for  a  5  micron  diameter  shell  for  several  thicknesses  are 
shown  in  Figure  1. 


Frequency  (GHz) 


Figure  1.  The  real  and  imaginary  parts  of 
the  permeability  as  a  function  of  frequency 
of  a  copper  shell  of  diameter  5  microns  with 
thicknesses  of  0.2S  microns  (Curves  1), 

0,10  microns  (Curves  2),  and  0.050 
microns  (Curves  3). 


It  was  hoped  that  some  insight  could  be  gained  from  the  calculated  behavior  of 
the  internal  field.  The  results  of  a  sample  calculation  for  the  real  and  imaginary 
parts  of  the  internal  field  for  a  5  micron  diameter  shell  0.25  microns-thick  are 
shown  in  Figure  2  (compare  with  Curves  2  in  Figure  1).  The  imaginary  parts  of 
the  internal  field  and  of  the  permeability  peak  at  essentially  the  same  frequency. 


Figure  2.  The  reel  end  imaginary  peris  of 
the  internal  fleW  as  a  function  of  frequency 
of  a  5  micron  diameter  copper  shell  of 
thickness  0.25  microns.  Similar  behavior 
is  obtained  for  other  shells. 
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COMPOSITE  BEHAVIOR 


Many  expressions  for  the  dielectric  and  permeability  constants  of  a  composite 
have  appeared  in  the  literature  over  the  years.  While  all  of  these  seem  to  have 
some  utility,  no  one  expression  seems  to  be  truly  adequate  for  all  circumstances.3'4 5 
One  of  the  most  commonly  used  expressions  is  that  of  Maxwell-Gamett,  which  is 
used  here  for  want  of  something  more  definitive.  Also,  the  interest  at  the  moment 
is  to  see  how  the  permeability  of  the  composite  depends  on  the  shell  thickness  and 
conductivity,  and  on  the  diameter  of  the  sphere  rather  than  on  the  volume  fraction 
of  the  spherical  particles.  In  this  approximation,  the  composite  permeability,  pt c,  can 
be  written  as: 


f*c  =  1 


+ 


3fpP-l) 

(l-O/p  +  2  +  f  ’ 


(21) 


where  f  is  the  volume  fraction  of  particles. 

Figures  3A  and  3B  show,  respectively,  the  real  and  imaginary  parts  of  the  perme¬ 
ability  as  a  function  of  frequency  of  a  composite  of  20  volume-percent  of  spherical 
particles  in  a  nonmagnetic  binder.  The  outer  radius  of  the  particles  was  2.5 
microns,  and  the  shell  thicknesses  are  as  indicated.  The  conductivity  of  the  shells 
was  taken  to  be  that  of  bulk  copper.  Similar  curves,  shifted  in  frequency,  are 
obtained  for  larger  diameter  shells  and  for  shells  of  different  conductivity.  For 
example,  results  for  the  real  and  imaginary  parts  of  the  permeability  of  50  micron 
diameter  copper  shell  are  shown  in  Figures  4A  and  4B,  in  which  a  log-linear  plot  is 
used  to  better  illustrate  the  behavior  of  the  imaginary  component. 

Two  interesting  features  appear  in  the  results  for  the  imaginary  component  of  the 
permeability.  First,  the  frequency  of  the  peak  in  ft”  (Fm„)  moves  up  in  frequency  as 
the  shell  thickness  decreases,  and,  secondly,  the  peak  value  of  pt”  also  increases  by  as 
much  as  35%  over  that  of  a  solid  particle.  A  plot  of  Fmax  as  a  function  of  shell  thick¬ 
ness  for  several  shell  outer  diameters  is  shown  in  Figure  5.  It  is  noted  that  there  is  a 
region  of  solid  particle-like  behavior  (in  which  Fmu  is  basically  independent  of  shell 
thickness)  and  of  shell-like  behavior  (in  which  F^  varies  inversely  with  shell  thick¬ 
ness),  and  that  the  transition  between  them  is  rather  abrupt  on  the  log-log  plot.  We 
also  find  that  for  a  given  shell  thickness,  FmM  varies  inversely  with  the  particle  diameter 
and  conductivity.  Finally,  for  very  high  frequencies,  a  =  -  3/8  n  and  =  0,  which  in 
turn  makes  a  pic  real  constant  less  than  unity. 

It  may  be  more  revealing  to  show  plots  of  rather  than  fic\  for  comparison,  a 
plot  of  Fm.»  for  /4p  is  included  in  Figure  5  for  a  5  micron  diameter  particle.  The 
behavior  is  similar  to  that  of  the  composite,  as  expected. 


4.  VANBEEK,  LK.H.  Dielectric  Behavior  of  Heterogeneous  Systems  in  Promt  in  Dielectrics.  J.  B.  Birkt,  cd., Chemical  Rubber 
Company  Press,  Cleveland,  OH,  v.  7, 1967,  p.  69-114. 

5.  HO,  Y.  S.,  and  KRAMER,  J.  J.  Microwave  Dielectric  Properties  of  Metal  Filled  Particulate  Composites  in  Microwave  Processing  of 
Materials.  W.  H.  Sutton,  H.  H.  Brooks,  and  I.  J.  Cha  bins  Icy,  ed.,  Materials  Research  Society,  Pittsburgh,  PA,  1988,  p.  161-166! 
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Rguree3Aand3B.  Th« r*al  and  Imaginary  parti  of  th* 
permeability  a*  •  function  of  frequency  of  a  oompoeit*  of 
20  vofcune-peroent  of  5  micron  diameter  oopper  aheHe  of 
thlokne****  2.5  (Le.,  a  eoUd  sphere),  0.5, 0.25, 0.1 , 0.06, 
and  0.025  micron*  (Curvee  1  -6) ,  re*p*o#v*iy.  The 
Maxweil-Oamett  axpraeaton  forth#  permeability  of  th# 
compoth*  w*a  u*#d  In  the  calculation*. 
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Figure*  4A  and  48.  The  real  and  imaginary  part*  of 
th*  permeability  a*  a  function  of  frequency  of  a 
eompoait*  of  20  volume-paroant  of  50  micron  diameter 
copper  ahaila  of  thickneea**  25  (aoiid),  5, 1,  and  0.1 
micron*  (Curve*  1-4),  raepactlvafy.  Notatha  log-Hnaar 
piot  Th#  Maxwaii-Qamatt  axpraaakxi  for  th*  permeability 
of  th*  compoaita  waa  uaad  in  th*  calculation*. 
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Figure  5.  A  plot  of  the  frequency  of  the  maximum 
in/i"  (Fmax)  veraue  ehell  thickneen  for  several  shell 
diameters  and  conductivities.  Tht>  Maxwell-Garnett 
expression  for  the  permeability  of  tni  composite 
was  used  in  the  calculations. 

Curve  1 :  S  micron  diameter,  copper  (particle) 
Curve  2:  5  micron  diameter,  copper  (composite) 
Curve  3:  10  micron  diameter,  copper  (oompoeite) 
Curve  4:  SO  micron  diameter,  copper  (oompoeite) 
Curve  5:  50  micron  diameter,  material  1/100  as 
conductive  as  copper  (composite) 


The  physical  reasons  for  most  of  these  results  are  not  readily  evident;  one  might 
have  expected  Fmax  to  peak  when  the  shell  thickness  is  roughly  equal  to  the  skin 
depth.  For  copper,  however,  the  skin  depth  at  10  GHz  is  0.66  microns,  whereas  Fmax 
is  about  10  GHz  for  shell  thicknesses  of  0.25,  0.1,  and  0.025  microns  for  5,  10,  and  50 
micron-diameter  shells,  respectively.  The  product  of  shell  diameter  and  Fmax  is  thus 
nearly  constant.  Since  the  shell  thickness  is  (generally)  appreciably  less  than  the  skin 
depth,  we  can  approximate  the  current  in  the  shell,  I,  as  being  uniform  within  the 
shell,  and  equal  to  2?rrtj,  where  r  is  the  radius,  and  t  the  thickness  of  the  shell. 

The  current  in  the  shell  and  the  internal  field  must  be  linearly  related,  and  from 
inspection  we  find  that  the  imaginary  part  of  the  internal  field  varies  linearly  with 
frequency  at  low  frequencies  (see  Figure  2).  Thus, 

fimag  =  2nTtrjaF ,  (22) 

where  rj  is  a  constant  and  the  conductivity  dependence  has  beer  made  explicit.  At 
very  high  frequencies,  the  imaginary  part  of  the  current  must  drop  to  zero  as  the 
penetration  depth  decreases,  and  Ijmag  will  therefore  peak  at  some  intermediate 
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value,  Fmax.  It  is  expected  that  the  maximum  in  Ijmag  ( iJiTtija  Fmax)  would  be  a 
constant,  from  which  the  dependence  of  Fmax  on  a,  r,  and  t  follows  directly. 

It  is  further  noted  that  expressions  for  a,  and  fic  are  mathematically  linear  in 
A  -  B,  and,  hence,  in  shell  thickness  for  small  shell  thicknesses. 

Other  models  for  the  dielectric  constant  of  a  mixture  of  spherical  particles  in 
a  matrix  could  also  be  extended  to  the  permeability.  For  comparison  purposes, 
Figures  6A  and  6B  show  the  expected  behavior  for  a  composite  of  20  volume 
percent  particles  calculated  from  the  Maxwell-Garnett  theory  (Curve  1),  from 
effective  medium  theory  (Curve  2),  for  which  is  given  by 


f  C“p  -f*c)  . 

+  2  nc 


(1-0 


(1  -Me) 
1  +  2  Me 


=  0  , 


(23) 


and  from  a  model  developed  independently  by  Looyenga  and  by  Landau  and  Lifschitz 
(Curve  3)  for  which 


[l+0*p1/3-l)]3  •  (24) 

The  real  part  of  the  permeability  (see  Figure  6A)  is  pretty  much  the  same  in 
the  models,  but  the  imaginary  part  (see  Figure  6B)  is  quite  different  at  high  frequen¬ 
cies. 


It  is  quite  possible  that  none  of  these  model  formulas  will  adequately  describe 
the  permeability  of  the  composite.  This  point  will  have  to  await  experimental  verifi¬ 
cation  (which  will  be  the  subject  of  a  subsequent  article).  However,  we  do  expect 
the  general  features  described  above  to  be  observed. 

The  dielectric  constant  of  a  composite  of  metal  particles,  or  shells,  in  an  insu¬ 
lating  binder  is  also  of  interest.  The  Maxwell-Garnett  expression  for  the  dielectric 
constant  of  such  a  composite,  ec,  simplifies  to 

e=  =  «b(  V=t0’  (25) 


where  ec  is  the  dielectric  constant  of  the  binder,  since  the  dielectric  constant  of  the 
metal  particle  is  extremely  large.  This  expression  seems  to  generally  underestimate  the 
dielectric  constant  and  such  microstructural  variables  as  particle  size  and  dispersion,  as 
well  as  material  parameters  (particle  type,  conductivity,  etc.),  are  ignored,  as  are  all 
effects  due  to  particle -particle  contact.  Ho  and  Kramer5  present  a  modified  expression 
for  the  dielectric  constant: 


[  ( 1  —  <h>)  (1  -  f)  +  <fa>  (I 
L  (1  —  £)  (1  —  <h>)  J  * 


(26) 


i 
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in  which  <h>  is  an  average  particulate  shape  factor,  which  serves  as  a  measure  of 
particle  contact  or  chain  formation,  which  is  obtained  by  an  independent  micros- 
tructural  analysis.  For  isolated  spherical  particles,  <h>  =  2/3,  and  approaches  1  for 
a  long  chain  of  contacting  particles.  Their  formula  gives  a  generally  improved  fit  to 
their  data  for  the  real  part  of  the  dielectric  constant  for  their  low  loss  samples. 


Figure*  6A  and  68.  The  real  (8A)  and  imaginary  (88)  part* 
cf  the  permeability  aa  a  function  of  frequency  ter  a  oompoelt* 
of  20  volume- par  cant  of  50  micron  diameter,  0.1  micron-thick 
spherical  sheila  calculated  using  different  models  for  the  permeability 
of  the  composite  (Curve  1 ,  Maxwall-Oamatt;  Curve  2,  Effective 
Medium  Approximation;  Curve  3,  Landau  and  Ufschitz  •  Looyenga). 
Mathematical  expreealons  for  the  permeability  are  given  in  the  text 
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Although  the  Ho-Kramer  article  does  represent  a  definite  improvement  in  the 
theory  of  the  electromagnetic  properties  of  metal-filled  particulate  composites,  it  is 
clear  that  more  improvement  is  desirable.  It  is  noted  that  Ho  and  Kramer  do  not 
present  data  for  the  permeability  of  their  samples;  this  is  largely  because  they  do 
not  understand  their  results  as  yet.* 


CONCLUSIONS 

An  expression  for  the  permeability  of  a  conducting  spherical  shell  which  should  be 
accurate  as  far  as  classical  electromagnetic  theory  is  concerned  has  been  derived.  We 
have  also  explored  the  calculated  behavior  of  composites  of  such  particles  in  an  insu¬ 
lating  binder  using  various  theories  for  the  permeability  of  the  composite.  A  subse¬ 
quent  article  will  address  a  comparison  of  the  calculated  results  with  experimental  data. 


•J.  J.  Kramer,  private  communication,  July,  1989. 
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